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■ Abstract. We consider some examples of quantum super- integrable systems and the as- 
sociated nonlinear extensions of Lie algebras. The intimate relationship between super- 

■ integrability and exact solvability is illustrated. Eigenfunctions are constructed through the 
^ , action of the commuting operators. Finite dimensional representations of the quadratic al- 
gebras are thus constructed in a way analogous to that of the highest weight representations 
of Lie algebras. 
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O ■ 1 Introduction 

Oh' In classical mechanics, there is a well defined meaning to the term 'complete integrability'. 

r~| ■ If, in n degrees of freedom, we have n independent functions in involution (mutually Poisson 

I commuting), then the system can be integrated, 'up to quadrature'. This is known as Liouville's 

Theorem. It is customary to consider one of these functions as 'the Hamiltonian' LI, with 
the others being its first integrals and this Hamiltonian is said to be completely integrable in 
the Liouville sense. Whilst n is the maximal number of independent functions which can be 
^ ' in involution, it is possible to have further integrals of the Hamiltonian LJ, which necessarily 

■ generate a non-Abelian algebra of integrals of H. The maximal number of additional independent 

integrals is n— 1, since the 'level surface' of 2n—l integrals (meaning the intersection of individual 
level surfaces) is just the (unparameterised) integral curve. Well known elementary examples 
are the isotropic harmonic oscillator, the Kepler system and the Calogero-Moser system. The 
quadratures of complete integrability are often achieved through the separation of variables of 
the Hamilton- J acobi equation. The solution of a maximally super-integrable system can also 
be calculated purely algebraically (albeit implicitly), requiring just the solution of the equations 
4 = Cfc, A; = 1, . . . ,2n - 1. 

In n-dimensions. Quantum integrable systems are defined analogously by requiring the exis- 
tence of n mutually commuting differential operators (usually with some requirements, such 
as self-adjointness). One of these will be of Schrodinger type, being the quantum version of 
the Hamiltonian H. However, we don't have a theorem analogous to Liouville's and we don't 
have anything resembling 'reduction to quadrature'. There is a notion of separation of variables 
(this time in the linear PDE sense), requiring a solution to be expressible as the product of n 
functions of a single variable. This reduces the problem of solving one Schrodinger equation in 
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n-dimensions to that of solving n Schrodinger equations in 1-dimension. This achieves less than 
quadratures. 

What we would like to achieve in quantum mechanics is to build the spectrum and corre- 
sponding eigenfunctions for Schrodinger's equation with some given boundary conditions. This 
class of Schrodinger equation is called exactly solvable. Even in 1-dimension exactly solvable 
Schrodinger equations are rare, so unlikely to be the end-product of applying the separation of 
variables method. One particular class of 1-dimension exactly solvable Schrodinger equations 
was presented by Infeld and Hull [8], related to the factorisation of operators. This can be 
related to Darboux transformations, which can then be generalised to higher dimensions (where 
factorisation is not possible). Darboux transformations can then be used to build exactly solv- 
able Schrodinger equations in higher dimensions. In [3], Darboux transformations have been 
shown to be related to the symmetries of the Laplacian. 

In the quantum case we also have a notion of super-integrability. By analogy with the classical 
case, super-integrability again involves having 'extra' commuting operators. Once again, we 
cannot have more than n such operators in involution, so the algebra is necessarily non-Abelian. 
Furthermore, if the operators are second order, their commutators are third order, with further 
commutators being fourth order and so on, so cannot be expected to generate a finite dimensional 
Lie algebra. However, they can generate a finite dimensional algebra with nonlinear commutation 
relations (see [211101 IT3]. for example). There are many papers on super-integrable systems (both 
classical and quantum) . Particularly relevant to the present paper is a series of papers by Kalnins 
et al. (see [H [10] and references therein), in which a classification is given of super-integrable 
systems in 2-dimensional flat and constant curvature spaces and in which a pair of quadratic first 
integrals, together with their quadratic commutation relations are derived. As is well known, 
super-integrable systems are separable in more than one coordinate system. In most of these 
papers the main emphasis is on the separation of variables, first by constructing and classifying 
coordinate systems and corresponding 'separable potentials' and then using this (in the quantum 
case) to build eigenfunctions. 

In the present paper, there is no discussion of separation of variables and 'differential equa- 
tions techniques' are kept to a minimum. Our main emphasis is on the direct relationship 
between super-integrability and exactly solvability. We use the commuting operators in a direct 
and explicit manner to build the eigenfunctions for our operators. This construction is the 
analogue (in the context of our quadratic algebras) to that of the highest weight representation 
of simple Lie algebras (see the Appendix). In particular, the polynomial eigenfunctions of the 
Krall-Sheffer operator of Example [9] can be explicitly constructed in this way. The issue of 
the relationship between super-integrability and exact solvability is not new and is discussed 
in [101 fill I14j . In [13] it is conjectured that exactly solvability should hold quite widely for 
super-integrable systems. Whilst operator algebras are discussed in these papers they are not 
used to construct eigenfunctions. 

In this paper we consider some particular super-integrable systems, associated with a par- 
ticular Laplace-Beltrami operator in 2-dimensions. One of the examples is in the Krall-Sheffer 
classification of 'admissible operators'. The relationship between Krall-Sheffer operators and 
super-integrability was shown by Harnad, et al in [6l [15]. In these papers they reworked the 
Krall-Sheffer classification and showed that in each case the leading order terms in the opera- 
tor correspond to the Laplace-Beltrami operator of either a zero or constant curvature, 2-di- 
mensional space. They also show that each of the 9 cases is super-integrable by presenting pairs 
of commuting, second order operators. Example [9] of the current paper is exactly Case II of [6lll5j. 
In [15] it is pointed out that this case had been previously introduced and studied in [Tl [TO]. 

In the next section, we give some basic formulae related to our particular metric and its 
symmetries. We then present 3 super-integrable examples, with their operator algebras and 
finite representation spaces (families of eigenfunctions). We then give some concluding remarks. 
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2 A metric and its symmetries 

For an n-dimensional (pseudo-)Riemannian space, with local coordinates and met- 

ric gij, the Laplace-Beltrami operator is defined by Lbf = g^-'Vi'Vjf, which has explicit form 



i,j=i vy 



where g is the determinant of the matrix g^j. The coefficients of leading order terms in the Lap- 
lace-Beltrami operator are the coefficients of the inverse metric g^K For a metric with isometries, 
the infinitesimal generators (Killing vectors) are just first order differential operators which 
commute with the Laplace-Beltrami operator ([1]). When the space is either flat or constant 
curvature, it possesses the maximal group of isometries, which is of dimension n(n + l)/2. In 
this case, is actually the second order Casimir function of the symmetry algebra (see [5]). 

In this paper, we consider one particular constant curvature metric in 2-dimensions, with 
inverse 

g'' = f V (2) 

\xy y^-yj 

The choice of coordinates here is motivated by the relationship to Krall-Sheffer operators |12] . 
but otherwise quite arbitrary. 

A convenient basis of Killing vectors is 

U = 4xd^, B = 2^dy, F = 4^d^ + 2{y-l)^dy, (3) 

satisfying the standard commutation relations of s/(2,C): 

[H,E]=2E, [H,F] = -2F, [E,F]=H. (4) 

The Laplace-Beltrami operator for the metric ([2D is proportional to the quadratic Casimir 
operator: 

L, = ^(H^ + 2EF + 2FE) = x^d^ + 2xyd,dy + {y^ - y)d^ + ^xd^ + ^(3y " l)^^- (5) 
Remark 1. The operator Lf^ is invariant under an involution, 

x = , y = y. (6) 

X 

which, in fact, gives a concrete realisation of the Lie algebra automorphism 

E ^ F, H ^ -H. 
This will be useful in some of the calculations below. 

Second order operators, commuting with Lf,, are just symmetric quadratic forms of Killing 
vectors. Suppose K is such an operator. Then we may seek functions u and v, such that 

[Lb + u,K + v] = 0, 

which constitutes a coupled system of partial differential equations for u and v. The solution 
depends upon a pair of arbitrary functions, each of one variable, typical of separable systems. 
Requiring that there exist two of these second order commuting operators, Ij = K.j+Vj, strongly 
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constrains these (formerly) arbitrary functions, which reduce to rational functions depending 
upon only a finite number of parameters. The coefficients of the second order derivatives in 
such an operator K define a contravariant, rank-two Killing tensor. For brevity, we will refer to 
these operators as Killing tensors in what follows. 

The examples discussed below were presented in [3], in the context of Darboux transforma- 
tions, which will not be discussed in this paper. Instead we discuss how to use the commuting 
second order operators to build eigenfunctions of L = -|- u in a way analogous to the highest 
weight representation of a Lie algebra. 

For some calculations, such as those of j3j, it is convenient to use separation coordinates. In 
our case we can choose 2 (or even 3) such coordinate systems, but the current calculations are 
as easily done in terms of the x — y coordinates. 

3 The system with KiUing tensors and 

Requiring both Ii = -|- vi and I2 = 'Ei'^ + V2 to commute with L leads to the specific forms of 
u{x,y), Vi{x,y) in the operators below 

L = Lb-\ 1 h C2 K—, 

X y x^ 

\ y X, 

(7) 

In order to build eigenfunctions of L it is convenient to change gauge: 

L^L = G-^LG with G = exp(^|(^^^^x('3+^-i)/V^27+i)/4^ (8) 

which simultaneously reduces the 3 potential terms to constants. The specific coefficients in G 
were chosen to simplify the final expressions for the operators. Removing the additive constants 
and adjusting overall multiplicative factors (and dropping the 'tilde') we arrive at the following 
operators: 

L = x^dl + 2xyd^dy + (y2 -y)dl + {fix + a)d.^ + {f3y + ^)dy, 
h = xydy + {ay - -/x)dy, 

l2 = x^dl + {{p + j)x + a{l-y))d^, (9) 




where cq = a{2 — (i — 7)/2, ci = (27 -|- 1)(27 -|- 3)/16, C2 = a^/4. This is an example of 
admissible operator, as defined by Krall and Sheffer [12]. The connection between Krall-Sheffer 
operators and super-integrability was first pointed out in [U [15] . Indeed, Example [9] above is 
exactly Case II of |6l[15]. As pointed out in [15], this case was previously introduced and studied 
in [llllO]. 

The operator = [/i, 12] is third order, so cannot be written as a polynomial (in particular, 
linear, Lie algebraic) expression in L, Ji, I2, but does satisfy the polynomial equation: 



H = ^{Iih + hhh + hit - 2It) + a(/3 - 7 - 2) 1^/1/2 + hh - 3/1 

+ (/3 + 7 - 2)((/? + 7)/? - 2ahL) + a\l2 - if + '^o?{l2 - (37 + l)L). 
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These 4 operators satisfy the commutation relations: 

[L,h] = [L,h] = [L,h]=0, [h,h]=l3, 
[hJs] = 2/f + a(/3 - 7 - 2)/i + a^h - L), 

[h, h] = -2{hh + hh) + (/? + 7 - 2)(aL - (/? + 7)^1) - - 7 - 2)/2. (10) 

By definition (see [E]), an admissible operator possesses a sequence of eigenvalues A„ (with 
An 7^ Am for m 7^ such that, for each n > 0, there exist n+1 linearly independent polynomial 
eigenfunctions of degree n, with eigenvalue A„. The degeneracy (for n > 0) of A^ stems from the 
super-integrability in a very explicit way. We show below how to use the operators Ii and I2 to 
build sequences of polynomial eigenfunctions of L, using a construction which is analogous to 
that of the 'highest weight representation' of a simple Lie algebra. 

We start with a polynomial eigenfunction 

n 

i=l 

for constants pi, which is independent of y and therefore in the Kernel of Ii. This satisfies the 
1-dimensional eigenvalue problem 

LxPnfl = {x^dl + {fix + a)d.j:)Pnfl = XnPnfi, A„ = n(n + /? - 1). 
This value of An is easily determined by looking at the coefficient of x". It is easy to show that 
a „ 2 2ax o? 



-Pi,o = 2; + -, ^2,0 = 2;^ + -^^ + 



/?' • /3 + 2 (/3 + 2)(/? + l)' 

and that this sequence satisfies the 3-point recursion relation 

'■ -fn+1,0 = {x + An)Pn,0 + BnPn-1,0, 

with 

_ a{P-2) _ a^n{n + (3-2) 



(/3 + 2n)(/3 + 2n-2)' (/3 + 2n - l)(/3 + 2n - 2)2(^ + 2n - 3) ' 

We now use I2 (for q 7^ 0) to build a sequence of [monic) eigenfunctions 

Pn~k,k = x"'~^y'' + lower order terms, k = 1, . . . ,n. 

The case of a = is the reduction cq = C2 = 0, which has a different algebra and is dealt with 
later (see the system (jl2p ). Since I2 commutes with L, all such eigenfunctions will have the same 
eigenvalue An, which means (important for Po.nj below) that the degree of each polynomial must 
be n. The polynomials Pn-k,k are defined recursively by 

l2Pn-k,k = {n-k){n-k-l+(3+-f)Pn-k,k-a{n-k)Pn-k-i,k+i, k = 0,...,n-l, (11) 

as can be seen by inspecting the form of I2 and by inspecting the coefficients of x^~^y^ and 
^n~k~iyk+i ^ It can be seen that /2^o,n = 0, so the polynomial Po,n must be independent of xJl 
The polynomials Po,n(y) also satisfy a 3-point recursion relation [3] 

: -Po,n+l = {y + -^n)PQ,n + BnPo,n-l-, 



^The condition (3.8) of [12], which ensures this, is that (3 should not be a negative integer 
^ Since the leading order term is killed, any non-zero result would be a polynomial of lower degree, which is 
not possible. This can also be seen by setting fc = n in (fTTj) 
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with 

(n + l)(7 — n) n(7 — n + 1) 



l3 + 2n /3 + 2ra-2 ' 



2(/3 + 2n-2) V '/3 + 2n-l /3 + 2n - 2 ' '/? + 2?i-3 
The space of polynomials spanned by this basis is also invariant under the action of /i: 
IiPn-k,k = kaPn-k,k + k{k--f - k = l,...,n. 

This array of polynomials, together with the action of Ii, I2 and the 3-point recursion relations, 
are depicted in Fig. [TJ 
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Figure 1. The triangular lattice of polynomials Pn-k,k, k = 0, ...,n with Po,n = 1, for system 
Horizontal arrows denote the action of Ii (left) and I2 (right). 

For each n > 1, the sequence of polynomials Pn-k,ki k = 0, . . . ,n, defines an n+1 dimensional 
representation of the algebra (jlOp . These can be explicitly calculated. For instance, apart 
from Pi^o and i-2,0, given above, we have 

7 1 Qf'Y 

-Po,i = y + -5, Pi,i =xy+ (jx + ay) + 



/?' " /3 + 2^' " (/J + 2)(/? + l)' 



P _„2 , n(7-l)y , 7(7-1) 
Po,2 — y + 2 ■ 



/3 + 2 (/3 + 2)(/3 + l)- 

It is easy to see that, for n > 2, the 4 representative matrices generate a larger Lie algebra 
(commutators lie outside the linear span), but do satisfy the quadratic relations ([TO]) . 

4 The system with Kilhng tensors and 

Requiring both Ii = + vi and I2 = F'^ + V2 to commute with L leads to the specific forms 
of u{x,y), Vi{x,y) in the operators below 

X ci x"^ 
L = Lt + Co 7 —5- H h C2 



/i = f2-4( ci^ — + C2 



'(y-l)2 y ^(y-1) 
{y - 1)2 



3 ■ 



(y-i)^, 

2^ 



/2 = H2 - 16 I CO ( ) + C2 



X 



Quantum Super-Integrable Systems as Exactly Solvable Models 



7 



These are, in fact, directly transformed from d?]) under the action of the involution Either 
by composing the gauge transformation ([S]) with this involution or by directly transforming Q 
by ([6]) (and dropping 'bars'), we obtain 

h = ^xydl + 4y(y - \)d^dy + ^^^^52 + 2 ((1 - 7)y + ^ + 7) 5, 
+ ay 5y, 



/2 = x2a2+((2-/3-7)x-^)9., 

where Co = a(2 — /? — 7)/2, ci = (27 + l)(27 + 3)/16, C2 = These satisfy the same quadratic 

relations (jlOp . Starting from these operators, it would not be so easy to find a simple function 
in the kernel of /i, and then to build the sequence of eigenfunctions, but there is no need to do 
this, since the functions 

/(y- 1)2 
Qn-k,k{x,y) = Pn-k,k ,y 



are obtained directly through the involution. The operator L is no longer 'admissible' in the 
Krall-Sheffer sense, and the eigenfunctions are no longer polynomial, but the same set of rela- 
tions hold, leading to the same diagram of Fig. [TJ 



5 The system with KiUing tensors and 

Requiring both /i = + vi and /2 = + V2 to commute with L leads to the specific forms of 
u{x,y), Vi{x,y) in the operators below 

L = L, + ^, /i = e2-4ci-, /2 = f2-4ci^^^. (12) 
y y xy 

However, it can be seen that the potential of L is a subcase of d?]), with cq = C2 = 0, so there 
should also be an integral of the form I3 = + ^3. Indeed, it should be exactly as in the case 
of d?!). However, in this reduction, we just have = H^, so may consider the algebra generated 
by L, Ii, I2 and H, which satisfy the commutation relations 

[H,/i]=4/i, [H,l2] = -4/2, [/i,/2] =2(8ci-1)H-H3 + 16HL, (13) 

as well as [L, Ii] = [L, Iq] = [L, H] = 0. Under the action of the involution ([6]), L ^ L, Ii ^ I2, 
I2 ^ /i, H — H, which preserves the commutation relations The Casimir operator of 

the algebra ^T3\i is 

C = hh + /2/1 + liSci - 5)H2 - Ih^ + 4U^L. 
2 

This is not, of course, independent of L and can be written C = 32L2-16(4ci + 1)L+16ci(2ci-1). 
We mimic the highest weight representation of sZ(2,C) and seek function ipi{x,y), satisfying 



Hyji = sipi, Lipi = Xipi, h<^i = 0. 



8 



A. P. Fordy 



The first of these gives ipi{x,y) = so the remaining equations lead to 

16y2(y - 1)^' + 8y{{s + 3)y - 1)'0' + (16ci + {s{s + 2) - 16X)y)tP = 0, 
2y2^' + yiP' - 2cii) = 0, 

which can be solved for both A and ij) to give 

Xm,s = T7^{s + m+l){s + in + ?,), ip 
lb 

We now define 

y^k = -f2~Vi, satisfying Ucp^ = {s - 4(/c - l))ipk. 
Under the involution, we have (pi i— > ipi and 

H(^i = l2<fi = 0. 

For this to be in our sequence tpk, we must have s — 4(A; — 1) = —s for some k. For some integer 
n, we therefore have 

s = 2n and (^i = i/3n+i, = 0. 

In this case, tpk, k = 1, . . . ,n + 1, define a finite dimensional representation of the algebra. For 
other values of s, the representation is infinite. 

The representations obtained in this way are deformations of representations of s/(2, C) (given 
in the Appendix), which correspond to the reduction ci = (which is m = — 1). In this reduction, 
the eigenvalue reduces to the value A_i^2n = \n{n + 1). When ci = 0, I2 reduces to F^, so our 
sequence ipi,(p2, ■ ■ ■ reduces to ip^jil^^,..., jumping in steps of 2. We thus retrieve only n + 1 of 
the 2n + 1 vectors of the highest weight representation. 

Example 1 (The case n = 2). Here we have 

ifi = xyP, ip2 = 4(3y - 1 + Ap{y - 1))/, ^3 = 8(y - 1)2(3 + 16p + Wp^)yP, 

where m = — 1 + Ap. When p = 0, these reduce to ipf, ijj^ and ip'^, given in the Appendix. 

6 Conclusions 

The main message of this paper is that super-integrability in a quantum system leads directly and 
explicitly to the construction of eigenfunctions. This leads to exact solvability, which is much 
stronger than just complete integrability. 

The examples presented here are fairly simple, but the general approach can be used for any 
super-integrable system. The examples presented here are associated with the Laplace-Beltrami 
operator of a 2-dimensional space of constant curvature. The main feature used was the existence 
of a large number of symmetries of the underlying space. Using these it is possible to deform the 
Laplace-Beltrami operator, by adding a potential, in such a way that it possesses (generally) 
higher order (non-geometric) symmetries in the form of higher order commuting operators. This 
can be thought of as a 'ghost' of the previous geometric symmetry. These higher symmetries 
generate an algebra with nonlinear commutation relations, whose representations are built in 
a way analogous to the highest weight representations of Lie algebras. In 2 dimensions the 
symmetry algebra is very small, being just s/(2,C) in our case (the 2-dimensional Euclidean 
algebra in the flat case), so the representations are very simple. In higher dimensions, the 
usual highest weight representations apply for the symmetry algebra of the Laplace-Beltrami 
operator. In this case there is no longer a single chain of eigenfunctions, so we can expect similar 
behaviour for the algebra of higher symmetries. 



y 



(m+i)/4^ where ci 



— 1 
16 
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Appendix. Eigenfunctions of the Laplace— Beltrami operator 

We use the highest weight representations of the symmetry algebra s/(2,C) to construct eigen- 
functions of Lb. Since H and Lb commute, they share eigenfunctions and for H these are built 
by the highest weight construction. Furthermore, starting with any eigenfunction of Lj, we may 
use the symmetry algebra to construct further eigenfunctions (with the same eigenvalue). Since 
this eigenspace is invariant under the action of sl{2, C) (by construction), it can be decomposed 
into irreducible components, which are just weight spaces. Therefore all eigenfunctions of L^ 
can be written as linear combinations of those we construct below. 
A highest weight vector ip'^, of weight 2n, satisfies 

which constitute a pair of partial differential equations for the eigenfunction. These are com- 
patible on the zeros of the differential operator E, since 

HEVi - EH^A]^ = 2E^^ = 0. 

The specific form of Vi depends upon n and upon the choice of representation for sl{2,C). 
However, the general structure of the representation is independent of this specific form, being 
a consequence only of the commutation relations (j3|) (see [7]). 
Defining ip^!: = F^~^ip^, the commutation relations imply: 

HV" = 2(n + l-r)V';?, E V'" = (r - l)(2?i + 2 - r) (14) 

Our definition of Laplace-Beltrami operator Lj, as Casimir operator ([5]) implies that 

Lbip^ = -^n{n + 1)^"^ , for all r, n. 

We can also construct a 3-point recursion relation between these eigenfunctions, but this is 
explicitly dependent upon the form of the operators. Let 

H = /ii(9^i +/i2 5^2' E = eia^, +62^22, F = fid^-, + f2d^^, 
where hi, etc are functions of the coordinates Zj. We have (with fir and defined by ([T 






hie2 - 61^2 
The relation V'r+i = F V'" then implies that 

in /ie2 - ei/2 „ /il/2 - /1/12 ,„ 
1pr+l = T 7-l^r1pr + T T-flrV'r-l- 

hie2 - ei/i2 hie2 - ei/12 

When this is singular the representation reduces to 1 dimension. 

In the case of our vector field representation ([3]) , this takes explicit form 

V-^i = 2(n - r + 1)./^C + {r- l)(2n + 2 - r) (^) 

with 



= 1^2 = 2n\/x'"-iy. 

When n is an integer, these representations are of finite dimension 2n + 1 and irreducible, 
but infinite dimensional otherwise. For instance, when n = 2, we have 

*| = 24fe-l),/^. *| = ?l<!i^. 
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